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Abstract. When two sheets of Color Glass Condensate collide in a high energy heavy 
ion collision, they form matter with very high energy densities called the Glasma. 
We describe how this matter is formed, its remarkable properties and its relevance 
for understanding thermalization of the Quark Gluon Plasma in heavy ion collisions. 
Long range rapidity correlations contained in the near side ridge measured in heavy 
ion collisions may allow one to directly infer the properties of the Glasma. 



1. Introduction 

At very high energies, multi-particle production in QCD is generated by low x partons 
in the nuclear wavef unctions. These partons have properties best described as a Color 
Glass Condensate (CGC) pQ. When two sheets of CGC collide in a high energy heavy ion 
collision, these partons are released and create energy densities an order of magnitude 
above the energy density required for the crossover from hadronic to partonic degrees of 
freedom. This matter, at early times after a heavy ion collision, is a coherent classical 
field, which expands, decays into nearly on shell partons and may eventually thermalize 
to form a Quark Gluon Plasma (QGP). Because it is formed by melting the frozen CGC 
degees of freedom, and because it is the non-equilibrium matter preceding the QGP, 
this matter is called the Glasma [2]. 

The Glasma is of intrinsic interest because heavy ion collisions create chromo- 
electric and magnetic fields in bulk, that in absolute magnitude, are some of the strongest 
such fields in nature; they are greater (by several orders of magnitude) than the magnetic 
fields on the surface of magnetars. The Glasma is also important for quantifying the 
initial conditions for the QGP. In particular, two measures of QGP formation, the flow of 
bulk matter and the energy loss of jets, can be significantly influenced by the properties 
of the Glasma. In the former case, the initial matter distribution and non-equilibrium 
flow are important for a quantitative determination of the properties of the "perfect" 
fluid QCD. In case of the latter, jets can experience significant energy loss in their 
interactions with the Glasma. 
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There is a strong analogy between the physics of the little bang in a heavy ion 
collision and the big bang that created our universe. In the big bang, the inflaton 
field with large occupation number 0{\) decays rapidly with the expansion of the 
universe. Likewise, in the little bang, the Glasma field with occupation number 
O(p-) decays rapidly after the collision. In the big bang, low momentum quantum 
fluctuations are explosively amplified in a process known as pre-heating [3J. In the 
little bang, the explosive amplification of low momentum quantum fluctuations may 
be related to a Weibel |] or a Nielsen-Olesen [5] type instability. The interaction of 
quantum fluctuations with the decaying inflaton field can lead to rapid "turbulent" 
thermalization [6]. A similar phenomenon may be responsible for rapid thermalization 
in heavy ion collisions [7] . There is also likely a concrete analogy between super horizon 
fluctuations observed in the COBE and WMAP measurements and the near side Ridge 
measured at RHIC [8J. Another strong analogy is between sphaleron driven topological 
transitions that may induce P and CP odd metastable states in heavy ion collisions [9] 
and the matter-anti-matter asymmetry generated by C and CP violating topological 
transitions during electroweak baryogenesis. What both long range rapidity correlations 
and topological transitions in strong fields have in common is the likelihood that they 
both survive the later stage interactions that thermalize the system causing it to lose 
memory of how it was formed. 

In this talk, I will first briefly describe our understanding of the properties of hadron 
and nuclear wavefunctions at high energies in the CGC framework. I will next describe 
the quantum field theory framework of particle production in strong time dependent 
fields. Multi-particle production in the Glasma can be computed systematically in this 
framework. An important part of this systematic computation is a proof of high energy 
factorization. I will describe briefly how plasma instabilities arise and are accounted 
for in this framework. Finally, I will discuss how Glasma flux tubes form the near side 
ridge seen in heavy ion collisions. 



2. Before the little bang 

At high energies, the competition between QCD bremsstrahlung which enhances the 
parton density)and screening/recombination processes which depletes it leads to a 
saturation of parton densities when the field strengths squared become maximal: 
E 2 ~ B 2 ~ This non-linear strong field regime of QCD is characterized by 

a saturation scale [H)] Qs{x,A); modes in the nuclear wavefunction with momenta 
k± < Qs have high occupation numbers typical of classical fields. In the CGC effective 
field theory, these classical fields are the dynamical degrees of freedom that couple 
stochastically to static light cone color sources p a (xj_) at large x [TT] : their source 
distribution is given by a gauge invariant weight functional Wy[p], where Y = ln(l/x) 
is the rapidity. While this separation of degrees of freedom is arbitrary at some initial 
scale, its evolution with energy is described by the renormalization group (RG) equation 
™^d£i. = H JIMWldK [p)W Y [p), where #jimwlk[p] is the JIMWLK Hamiltonian [12]. The 
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Balitsky-Kovchegov (BK) equation [T3] is a useful mean field (large N c , large A) 
simplification of the JIMWLK equation describing, in closed form, the energy evolution 
of the "dipole" operator corresponding to the forward scattering amplitude in deep 
inelastic scattering. The JIMWLK and BK equations are derived in the leading 
logarithmic approximation in x, where running coupling effects are neglected. There 
is currently an intense on-going theoretical effort to compute next-to-leading order 
corrections to the kernels of the leading order RG equations [14]. 

Phenomenological "dipole" models that incorporate the physics of saturation have 
been very successful-see Ref. [15] for a recent review of comparisons of these models 
to the HERA, fixed target e+A, D+A and A+A RHIC data. A detailed dipole model 
study [16] shows that the saturation scale Q 2 s (x,A) = A l f z Q\(x, A = 1). See fig. [I] for 
partial results of this study. Because of the nuclear "oomph" , saturation (CGC) effects 
should become increasingly visible in deuteron+gold collisions at RHIC, p+A collisions 
at the LHC, and especially at a future electron ion collider (EIC) |17j . 




Figure 1. Left: Parameters fit to HERA inclusive data in the CGC motivated IPsat 
and bCGC models give excellent agreement with HERA diffractive and exclusive vector 
meson data. As shown, the same models, with no additional parameters, compared to 
fixed target NMC data; rightmost figure contains curves for fixed Q 2 = 5 GeV 2 . See 
Ref. [16] for further details. Right: The saturation scale for a quark probe (for gluons 
multiply Qs by color factor |) as a function of x and A. The nuclear saturation 
scale extracted agrees with the value extracted from A+A multiplicity [20] at RHIC 
to - 10%. 

3. Multiparticle production in the Glasma at leading order 

When two CGC sheets collide in an A+A collision, the static color sources in the 
nuclear wavefunctions become time dependent. An ab initio computation of the 
Glasma therefore requires a formalism to compute particle production in Quantum Field 
Theories coupled to external time dependent sources [HI [19]. Multi-particle production 
in the Glasma is always non-perturbative; the relevant question is whether the physics 
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is one of strong coupling or weak coupling as will be considered here^l 

For the inclusive gluon distribution, the leading order (LO) contribution is of order 
0(^2 ) but all orders in gp\$- It can be expressed as 

E pyr = ttti lim I d3xd3 y elHx ~ y) ( d ° - iE p) ( d °y + iE p) 

xi^ZipKip) (M*)Mv)) ■ (i) 

A 

The gauge fields on the right hand side can be computed numerically for proper times 
t > [21] by solving the classical Yang-Mills (CYM) equations in the presence of the 
light cone current J M,a = 5^ + 5(x~)pl(x±) + 5^~5(x + )p%(x±) corresponding to the local 
color charge densities of the two nuclei with initial conditions determined by matching 
the (known) solutions to the CYM equations in the backward light cone. The energy 
densities of produced gluons can be computed in terms of Qs and one obtains e ~ 20 - 
40 GeV/fm 3 for the previously mentioned values of Qs obtained by extrapolating from 
fits to the HERA and fixed target e+A data. 

This LO formalism was applied to successfully predict the RHIC multiplicity at 
y ~ [21] as well as the rapidity and centrality distribution of the multiplicities [22J. 
CGC model comparisons to the RHIC data on limiting fragmentation [23] or solutions 
of CYM equations [20], extrapolated to the LHC, give dN/dri\ v ^ « 1000 - 1400 charged 
particles^. At LO, the initial transverse energy is E? ~ Qsi which is about 3 times larger 
than the final measured Et, while (assuming parton hadron duality) iVcGC ~ A^ad.- The 
two conditions are consistent if one assumes nearly isentropic flow which reduces Et due 
to PdV work while conserving entropy. This assumption has been implemented directly 
in hydrodynamic simulations |25j. 

CGC based models give values for the initial eccentricity e that are large than those 
in Glauber models [26J because the energy and number density locally is sensitive to the 
lower of the two saturation scales (or local participant density) in the former and the 
average of the two in the latter. Naively, CGC initial conditions would have more flow 
then and have more room for dissipative effects relative to Glauber. This conclusion is 
turned on its head in a simple parametrization of incompletely thermalized flow [27J: 

V2 l e = (i+ jr/ifo) ' wnere K = 'iz a ^ Cs * s ^ ne Knudsen number, a the cross-section, c s 
the sound speed and S± the transverse overlap area. Kq is a number of order unity. If 
thermalization were complete, K — > and one approaches the hydro bound. Computing 
e with different initial conditions, and plotting the l.h.s ratio versus s^^r, one has a two 
parameter fit to a and c s . The greater CGC eccentricity forces v^/eto be lower for more 
central collisions thereby leading to lower c s ; quicker saturation of v%/e forces larger a 
and therefore lower rj in the CGC relative to Glauber [28]. It is conceivable however 
that this result may not prove robust against more detailed modeling. Nevertheless, it 
is clear that the results are very sensitive to the initial conditions. 

| The magnitude of the saturation scale is Q 2 S ~ 1 - 1.4 GeV 2 at RHIC; estimates for the LHC are 

Q% ~ 2.6 - 3.9 GeV 2 [20]. 

§ See Ref. [24] for other model predictions. 
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How much flow is generated in the Glasma before thermalization? The primordial 
Glasma has occupation numbers f ~ — and can be described as a classical field. As 
the Glasma expands, higher momentum modes increasingly become particle like and 
eventually the modes have occupation numbers / < 1, which may be described by 
a thermal spectrum. A first computation of elliptic flow of the Glasma found only 
about half the observed elliptic flow [29] albeit the computation did not properly treat 
the interaction between hard and soft modes in the Glasma. Formulating a kinetic 
theory that describes this evolution is a challenging problem in heavy ion collisions-for 
a preliminary discussion, see [30J. 

The LO Glasma result, from the solutions of Yang-Mills equations, has very 
interesting properties. Firstly, the solution is boost invariant in the strong sense-the 
fields are independent of the space-time rapidity; the dynamics of the produced gluon 
fields is purely transverse as a function of proper time. Another interesting feature is 
that the chromo E and B fields are purely longitudinal after the collision [2]. This 
result suggests that one can generate topological Chern-Simons charge in heavy ion 
collisions [31]. Because the range of color correlations in the transverse plane is of order 
1/Qs, the LO picture that emerges is of color flux tubes with finite topological charge, 
stretching between the valence color degrees of freedom. As we shall see, this picture 
provides a plausible explanation of the near side ridge. 




Figure 2. Left: Chern-Simons mean squared charge as a function of proper time (in 
units of the saturation scale), generated in the LO boost-invariant 2+1-D Glasma [31] , 
Quantum fluctuations allow sphaleron transitions which may give rise to significantly 
greater values of yj {v 2 ). Right: The same quantum fluctuations, albeit suppressed by 
as, grow as as exp(y^gr) |32j . These unstable fluctuations, resummed to all orders, 
likely give rise to more isotropic initial distributions in the Glasma [19] . 

The LO field configurations are very unstable and lead to very anisotropic 
momentum distributions at later times. Such distributions can trigger an instability 
analogous to the Weibel instability in QED plasmas [I]. Romatschke and I showed 
(in 3+1-D numerical solutions of CYM equations) [32] that small rapidity dependent 
quantum fluctuations in the initial conditions generate transverse E and B fields that 
grow rapidly as ey^WQ s r ) ■ They are the same size as the rapidly diluting longitudinal 
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E and B fields on time scales of order -^-ln 2 (— ). The transverse E and B fields 

Qs v «s' 

may cause large angle deflections of colored particles leading to px broadening and 
energy loss of jets-numerical simulations by the Frankfurt group appear to confirm this 
picture [33J. These interactions of colored high momentum particle like modes with 
the soft coherent classical field modes may also generate a small "anomalous viscosity" 
whose effects on transport in the Glasma may mask a larger kinetic viscosity [31] . The 
same underlying physics may cause "turbulent isotropization" by rapidly transferring 
momentum from soft "infrared" longitudinal modes to ultraviolet modes \7\. Finally, 
albeit the LO result demonstrated that one could have non-trivial Chern-Simons charge 
in heavy ion collisions, the boost invariance of CYM equations disallows sphaleron 
transitions that permit large changes in the Chern-Simons number [31J. With rapidity 
dependent quantum fluctuations, sphaleron transitions can go. These may have 
important consequences-in particular P and CP odd metastable transitions that cause 
a novel "Chiral Magnetic Effect" §3] in heavy ion collisions. Numerical CYM results for 
Chern-Simons charge and (square root) exponential growth of instabilities are shown in 

fig.H 

4. QCD Factorization and the Glasma instability 

The discussion at the end of the last section strongly suggests that next-to-leading 
order (NLO) quantum fluctuations in the Glasma, while (superficially) parametrically 
suppressed, may alter our understanding of heavy ion collisions in a fundamental way. 
To cosmologists, this will not come as a surprise-quantum fluctuations play a central 
role there as well. In recent papers, it was shown for a scalar theory that moments of the 
multiplicity distribution at NLO in A+A collisions could be computed as an initial value 
problem with retarded boundary conditions; this framework has now been extended to 
QCD [18]. In QCD, the problem is subtle because some quantum fluctuations occur in 
the nuclear wavefunctions and are responsible for how the wavefunctions evolve with 
energy; others contribute to particle production at NLO. Fig. [3] illustrates particle 
production in fields theories with strong sources and the non-factorizable quantum 
fluctuations that are suppressed in the leading log framework. 

A factorization theorem organizing these quantum fluctuations shows that all order 
leading logarithmic contributions to an inclusive gluon operator O in the Glasma 
gives [19] 

(OK og = J [Dp,] [Dp 2 ] W Yh _. Y [Pi] W Yheam+Y [M O h0 [p h p 2 ] , (2) 

where O hQ is the same operator evaluated at LO by solving classical Yang-Mills 
equations and WY beam ^.Y [px,2\ are the weight functionals (introduced in section [2]) that 
obey the JIMWLK Hamiltonians describing the rapidity evolution of the projectile 
and target wavefunctions respectively. This theorem is valid if the rapidity interval 
corresponding to the production of the final state, AY < — . The Ws are 
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analogous to the parton distribution functions in collinear factorization; determined non- 
perturbatively at some initial scale Y , their evolution with Y is given by the JIMWLK 
Hamiltonian. 

This factorization theorem is a necessary first step before a full NLO computation 
of gluon production in the Glasma. Eq. ([2]) includes only the NLO terms that are 
enhanced by a large logarithm of \jx\ i, while the complete NLO calculation would also 
include non enhanced terms. These would be of the same order in ag as the production 
of quark- ant iquark pairs [35] from the classical field. To be really useful, this complete 
NLO calculation likely has to be promoted to a Next-to-Leading Log (NLL) result by 
resumming all the terms in asips ln(l/a;i j 2) ri . Now that evolution equations in the 
dense regime are becoming available at NLO, work in this direction is a promising 
prospect [T4] . 



In addressing the role of instabilities at NLO, note that small field fluctuations fall 
into three classes: i) Zero modes (p v = 0) that generate the leading logs resummed 
in eq. [2j the coefficients of the leading logs do not depend on x^. ii) Zero modes 
that do not contribute at leading log because they are less singular than the leading log 
contributions. These become relevant in resumming NLL corrections to the factorization 
result [Hj. Because they are zero modes, they do not trigger plasma instabilities, iii) 
Non zero modes (p v ^ 0) that do not contribute large logarithms of l/xxp, but grow 
exponentially as exp(y/Q s r). While these boost non-invariant terms are suppressed by 
as, they are enhanced by exponentials of the proper time after the collision. These 
leading temporal divergences can be resummed and the expression for inclusive gluon 
operators in the Glasma revised to read 



where Af(x) = — -^rPi(xj_,x ) and A 2 (x) = —-^rP2(x±,x + ). The effect of the 
resummation of instabilities is therefore to add fluctuations to the initial conditions of 




Figure 3. Left: Illustration of particle production in a field theory with strong time 
dependent sources. Right: Quantum fluctuations, where the nuclei talk to each other 
before the collision, are suppressed; this is responsible for high energy factorization of 
inclusive gluon operators O in A+A collisions [TO] , 




(3) 
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the classical field, with a distribution that depends on the outcome of the resummation. 
This spectrum Z[a(u)] is the final incomplete step in determining all the leading singular 
contributions to particle production in the Glasma-see however Ref . [36] . The stress- 
energy tensor can then be determined ah initio and matched smoothly to kinetic 
theory or hydrodynamics at late times. 

5. Two particle correlations in the Glasma and the Ridge 

Striking "ridge" events were revealed in studies of the near side spectrum of correlated 
pairs of hadrons at RHIC [8]. The spectrum of correlated pairs on the near side of the 
STAR detector extends across the entire detector acceptance in pseudo-rapidity of order 
Ar] ~ 2 units but is strongly collimated for azimuthal angles A<p. Preliminary analyses 
of measurements by the PHENIX and PHOBOS collaborations corroborate the STAR 
results. In the latter case, the ridge is observed to span the PHOBOS acceptance in 
pseudo-rapidity of Ar] ~ 3 — 4 units. 




Figure 4. Lcft:Glasma flux tubes of transverse size J- < -r-^ — ■ The field lines 

° Qs Aqcd 

correspond to parallel E and B fields, which carry topological charge. Right: The 
leading two particle contribution. Superficially disconnected, they are connected by 
averaging over the large x color sources. Systematic power counting shows these graphs 
dominate over usual "pQCD" graphs at high energies. 

Causality dictates (in strong analogy to CMB superhorizon fluctuations) that 
long range rapidity correlations causing the ridge must have occured at proper times 
t < Tf rccze out e~ 2 \va ~Vb I ) where y A and y B are the rapidities of the correlated particles. 
If the ridge span in psuedo-rapidity is large, these correlations must have originated 
in the Glasma. As noted previously in section [3J particles produced from Glasma flux 
tubes are boost invariant. See fig. HI Correlated two gluon production in the Glasma 
flux tube can be shown to be independent of rapidity [37]. Ours is the only dynamical 
model with this feature-for other models, see Ref.[38j. The particles produced in a flux 
tube are isotropic locally in the rest frame but are collimated in azimuthal angle when 
boosted by transverse flow [39]. Combining our dynamical calculation of two particle 
correlations with a simple "blast wave" model of transverse flow, we obtain reasonable 
agreement with 200 GeV STAR data on the amplitude of the correlated two particle 
spectrum relative to the number of binary collisions per participant pair. This simple 
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model also gives reasonable agreement for the amplitude for collisions at 62 GeV [40J. 
One caveat is that the collimation in azimuthal angle is weaker than seen in the data- 
this can be significantly improved with a more refined treatment of transverse flow. 
The Glasma flux tube model has several additional attractive features consistent with 
observations-see Ref. [37] for a fuller discussion. 
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